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1 Introduction 



The Lame equation, 

[E-mi{i + l) sn^x]ip{x) = 0, (1.1) 
where £ is a real parameter^, and 

sn X = sn(a;|m) 

is the usual Jacobian elliptic function of modulus m, occupies a central po- 
sition in the theory of differential equations with periodic coefficients. The 
study of its properties has attracted the attention of many illustrious mathe- 
maticians over the last century; classical references are |1|^ • Basic properties 
of the Lame equation are as follows. First, it arises by separation of variables 
in the Laplace equation in ellipsoidal coordinates. Secondly, it possesses two 
linearly independent 2K{k) or 4i^(A;)-periodic solutions (for characteristic 
values of E) if and only if £ is a nonnegative integer. Here k = \pm^ and 
K{}z) (denoted by K from now on) is the complete elliptic integral of the 
first kind with parameter k: 

K{k) = r'^ . 

Jo V 1 — /c^ sin^ X 

Moreover, if £ is a nonnegative integer there are exactly i gaps in the energy 
spectrum, and the 2i + 1 eigenfunctions associated to the boundaries of the 
allowed energy bands can be written as homogeneous polynomials of degree i 
in the elliptic functions sn, cn, dn. These polynomial solutions are known as 
the Lame polynomials. In the third place, for any £ G M and characteristic 
values of E the Lame equation admits two linearly independent solutions of 
period 8K. These solutions are expressible in closed form if and only if 
i G N + 1/2, and shall be referred to as the non-meromorphic Lame 
functions. For both types of solutions, namely the Lame polynomials for 
integer i and the non-meromorphic Lame functions for half-integer i, the 
characteristic values of E are the solutions of a certain algebraic equation. 
We shall use the the term "algebraic" to refer to both types of solutions. 

The interest of the Lame equation is by no means restricted to its mathe- 
matical properties. In fact, several important physical applications have been 
recently proposed in the literature. As first remarked in , the Lame poten- 
tial can be considered as a realistic model of a one- dimensional crystal. The 
Lame equation also describes some class of fluctuations of the sphalerons in 

-^Note that without loss of generahty we can assume that £ > —1/2. 
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the Abelian Higgs model in 1+1 dimensions, The quantum fluctuations 
of the inflation field in certain cosmological models are determined by the 
Lame equation, 

The Lame equation appears in a natural way in several Lie-algebraic 
approaches to the Schrodinger equation. Early work in this direction was 
carried out by Alhassid, Giirsey and lachello in 0. In this paper the Lame 
equation was obtained by separation of variables in the eigenvalue equation 
for a Hamiltonian quadratic in the generators of su(2) written in conical co- 
ordinates. A different Lie-algebraic representation of the Lame equation was 
discovered by Turbiner in where it was shown that for integer C. the Lame 
Hamiltonian belongs to the the enveloping algebra of a Lie algebra of first- 
order differential operators with a finite-dimensional module of functions. It 
follows that a number of eigenvalues and the corresponding eigenf unctions 
of the Lame equation can be determined algebraically, i.e., by diagonalizing 
the finite-dimensional matrix representing the action of the Hamiltonian in 
the module. This class of Hamiltonians — known in the literature as quasi- 
exactly solvable (QES) — have been the subject of considerable investigation 
over the last decade; extensive reviews of this field can be found in [|T0| -[12 



In particular, it is known that several periodic potentials of physical inter- 
est — like the Razavy potential, |13[|1^, or the so-called associated Lame 
potential, |Q — possess the QES property^. 

In the above mentioned papers the Lame equation was successfully alge- 
braized only for integer values of £, corresponding to the Lame polynomials. 
A first indirect algebraization of the Lame equation for half-integer i was 
studied by Ward in ||15[ in connection with the matrix Nahm equations. The 
Lame polynomials and the non-meromorphic Lame functions appear as solu- 
tions of some matrix- valued differential equations related to irreducible repre- 
sentations of so (3). The Lame potential (and some generalizations thereof) 
was revisited by Ulyanov and Zaslavskii in DTG). In this paper the Lame 
Hamiltonian was obtained from a non-standard realization of su(2) by first- 
order differential operators, leading to the algebraic solutions for integer and 
half-integer values of i. In [[T^, Brihaye and Godart introduced yet another 
indirect algebraization scheme for half-integer values of i based on a system 
of two second-order differential equations. 

In this paper we shall present a new direct way to algebraize the Lame 
equation which is valid for both the integer and half-integer cases. Our 
approach is based on the classification of one-dimensional QES potentials 
by Gonzalez- Lopez, Kamran, and Olver, |jTO|, and presents some advantages 
over the similar algebraization considered by Ulyanov and Zaslavskii in 

notorious exception is the Mathicu potential, which is not QES. 
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Indeed, the 2i + 1 algebraic points of the spectrum (corresponding to the 
boundaries of the energy bands for integer i and to eigenfunctions of quasi- 
momentum 7t/{4K) for half-integer i) are the roots of a pair of polynomials 
of degrees i and i + 1 for integer i, or of a single polynomial of degree 
£ + 1/2 for half-integer i. This compares to the algebraization by Ulyanov 
and Zaslavskii, for which one needs to find the roots of a polynomial of degree 
2£+l. In the second place, it leads to new general expressions for the algebraic 
solutions of the Lame equation in terms of Chebyshev polynomials and of a 
certain family of weakly-orthogonal polynomials. It should be emphasized 
that in the algebraization schemes mentioned above explicit solutions are 
presented only for low (integer or half-integer) values of i. 

The paper is organized as follows. In Section 2, we introduce the ba- 
sic definitions and show how the Lame potential can be obtained using the 
classification of Gonzalez-Lopez, Kamran, and Olver. The algebraization 
by Ulyanov and Zaslavskii is briefly discussed in this setting. Section 3 is 
devoted to the integer case. We find an expansion in terms of Chebyshev 
polynomials which reveals the structure of the classical Lame polynomials. 
This expansion is used to compute the Lame polynomials for i < 4. In the 
last section we analyze the half-integer case. We obtain an expansion of the 
non-meromorphic Lame functions using again Chebyshev polynomials, and 
relate our formulae to the classical results of Ince. In particular, we check 
Ince's formulas for i < 3/2, and study in detail the case i = 5/2. 

2 Algebraizations of the Lame equation 

If n is a non-negative integer, it is a standard fact that the differential oper- 
ators 



are the basis of a representation of the Lie algebra sl{2, M) in the space Vn of 
polynomials of degree at most n in the variable z, [|1^,|19|. As a consequence, 
any differential operator 7i which is a polynomial in the operators (|2.1|) pre- 
serves the space Vn (the converse is also true; cf. [^^). In particular, n + 1 
eigenfunctions and eigenvalues of H can be computed algebraically, simply 
by diagonalizing the finite-dimensional operator obtained by restricting 7i to 
Vn- More generally, if we perform the change of variable 



X 




J+ = z'^dz — nz 



(2.1) 



z = C{x) 



(2.2) 
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and consider the operator 



(2.3) 

z=ax) 



where fi{z) is a nowhere vanishing function, then H obviously leaves invariant 
the (n + 1) -dimensional vector space 

Consequently, we can again compute n + 1 eigenvalues and eigenfunctions of 
H in an algebraic way, by considering its restriction to the finite-dimensional 
space nVn- We shall say that such an operator H is quasi-exact ly solvable 



(or, for short, QES), |2^. Of particular physical interest is the case in which 
the operator 7i is a polynomial of degree two in the generators which 
we shall write as 

= - ^ Cab JaJb - ^ CaJa " . (2.4) 
a,fe=0,± a=0,± 

Such an operator can be expressed as 

-n = P{z) dl + [Q{z) -\{n-l) P'{z)\ d^ + R- -Q\z) + ^(n - l)P'\z) , 

(2.5) 

where (see |T9[] ) 

P{z) = C++ z^ + 2c+o z^ + (2c+_ + Coo) -^^ + 2co_z + c_._ , 

Q{z) = c+z^ + cqz + c_ , (2.6) 

R = Y^{n + 2) Coo - 3 ('^ + 2) c+_ + c, . 
Note that, using the Casimir identity 

we can take c^ = without loss of generality. It can be shown, |T^, that 

if P is positive then it is always possible to transform (at least locally) the 
operator 7i into a Schrodinger operator (or Hamiltonian) 

H = -dl + V{x) (2.7) 
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by an appropriate change of variable and gauge transformation ( |2.2| )- 
(|2.3| ). More precisely, the inverse of the change of variable ( p.2|) is given 

by 



X 



dy 



while the gauge factor is proportional to 



{2.1 



exp 



Q{y) 
my) 



dy 



(2.9) 



The potential V{x) can be expressed in terms of the coefficients of the gauge 
Hamiltonian Ti as follows 



V{x) 



12 P 



[n{n + 2) [PP" - IP'^) 

+ 3(n + 1)(QP' - 2PQ') - 3Q^] L^^^^^ - i?, (2.10) 



where the prime denotes derivative with respect to z. 

All one-dimensional potentials whose Hamiltonian is QES (in the precise 
sense explained above) have been completely classified in Ref. [|^, modulo 
a constant translation of the space coordinate x. In particular, the Lame 
equation ( |1.1| ) can be obtained from the third family of periodic potentials 
(with = 1), given by 



V{x) = A sn^ a; + i? sn a; cn a; + 



C sn X cn a; + D 
dn^ X 



[2.11) 



Note that, for later convenience, the above notation is slightly different from 
that of Ref. . The potential ( p.ll|) is obtained from ( |2.10|) when 



Piz) = (l + z^)[l + {l-m)z^]. 



Indeed, from ( |2.8| ) it follows that, in this case, 

sn X 



cnx 



(2.12) 



(2.13) 



Substituting (|2l^) into (plOD and using ( p.l3|) we find the following explicit 
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formula for the coefficients in Eq. ( p.llj ) 



A = -mnin + 2) 
4 ^ ' 



B = — (c+ - c_) [m{n + 1) - Cq] 
2m '- 



m 



l)c_] [m(n + 1) + Co] 



D = -(m - l)n(n + 2) + — (m - l)(n + 1) 



+ 



4m^ 



(m - l)co + (c+ + (m - l)c_)' 



(2.14) 
(2.15) 
(2.16) 

(2.17) 



Note that we have taken 
1 



4m^ 



[(2m -l)cl + (1 



+ 2c+c_) - 4] 



1 



n{n + 2) (2.18) 



in order to ehminate a constant term in the potential. Comparing ( |1 . Ij ) with 
(|2.11| ) we obtain the system 

A = me{e+1), B = C = D = 0, 

which has the following four sets of solutions: 



n 
n 

n 
n 



1; 
1 

2' 
1 



c+ 



c_ = 0, Co = — m£ 

c_ = 0, Co = -m (^ + 1) 



m, 



Co 



c_ 



-iv^ 



Co 



-m (£ + i) 
+ 1) 



(2.19) 
(2.20) 

(2.21) 
(2.22) 



The first two solutions are valid when £ is a non-negative integer {i > 1 
for the second solution to exist), while the last two solutions hold when i 
is a positive half-integer. From the previous general discussion it follows 
that when £ is a non-negative integer or a positive half-integer a certain 
number of eigenf unctions and eigenvalues of the Lame Hamiltonian can be 
computed in a purely algebraic fashion, by solving the eigenvalue problem 
of the restriction of the operator Ti. corresponding to the Hamiltonian ( |1 . 1| ) 
to the finite-dimensional space Vn, where n is given by ( |2.19|) - (|2.22|) . An 
explicit expression for the gauge Hamiltonian Ti. can be easily found from 
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Eqs. (|23), (gD, (Q, and {^)-(^), namely 



-n 



m 



m 



4 + (2 - m) 4 



Ji 



(1 - m) J2 



(2 
(2 



m) 
m) 



J_ — m £ Jo + c* 

J2 - m (£ + 1) Jo + c* 



J2 +iv/r^^ (J+ + J-) 

- m {i + ^) Jo + c* 
(1 - m) J2 + (2 - m) J2 + J! - i yr^^ ( J+ + J-) 

+1) Jo + c. 



(2.23) 



m 



Each of the four cases in this formula corresponds to the respective solu- 
tion ( |2.19| )- (|2.22|) , and it is understood that in each case the operators J^ 
(e = 0, +) and the constant c^: must be computed using the value of n for 
the corresponding solution ( p.l9| )- p.22| ). Note that Eq. (|2.23|) is essentially 
equivalent to Eqs. (7)-(9) of Ref. P|. The latter equations, however, ap- 
pear more complicated than ( |2.23| ), due to the fact that Ref. |0 uses the 
Weierstrassian form of Lame's equation. 

Before studying in more detail the four algebraizations ( |2.19|) - (|2.22| ) of 
the Lame equation obtained in this section, it is worth pointing out that there 
are other alternative algebraizations within the general formalism described 
above. Indeed, from the identity 



dn X 



m 



1 — m 1 — m 



sn^ X . 



m 



and the fact that the potentials in Ref. [|T0[ are classified up to constant 
translations, it can be shown that the Lame potential is a member of the first 
two families of periodic QES potentials listed in the latter reference. For the 
same reason, we could also have derived the Lame equation by equating to 
zero the coefficients A, B and C in Eq. ( |2.11| ). It can be shown that all these 
ways of algebraizing the Lame equation are essentially equivalent to the one 
adopted in this paper. Finally, there are other less obvious algebraizations 
of the Lame equation like, for instance, the non-standard one discussed in 
Ref. [ITH], which falls into our framework by taking 



n 



m 



J-] 



From Eq. ( |2.6| ) it follows that in this case 



P(z) = - 



m 



1 



It is straightforward to show that the complex change of variable 



n{n + 2) 
12 



cn x -|- i sn a; 
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(cf. Eq. (|2.8|)) and the gauge transformation defined by ( pl3| ) and ( |2.9| ) map 
7i into the Schrodinger operator (|2.7|) with potential 



T^/ \ .n/n \ sn^ a; 

V{x) = -m{l-m)-i^- + lJ 



dn^ a; 



= m|(^ + l) [sn2(x + K)-l] , (2.24) 

which is essentially the Lame potential with i = n/2. Since P{z) is negative 
everywhere and all its roots are complex, it can be mapped to a negative 
multiple of the polynomial ( p.l2| ) by a projective change of variable and 



gauge transformation. As explained in Ref. |T0|, this means that the potential 
(|2.24| ) can be obtained by a suitable choice of parameters in the third family 
of QES periodic potentials (with pure imaginary frequency y^)- 

3 Case I: ^ is a non-negative integer 

We shall show in this section that when £ is a non-negative integer the alge- 
braic eigenfunctions of the Lame Hamiltonian are the classical Lame poly- 
nomials, 0]. 

From Eqs. (|J), ( ^ ), and (|2l9|)- (|2:20|) we deduce that in this case the 
gauge factor verifies 

f^(cix)) = hl':; ^ ""=1 ^ (3.1) 

cn a; dn X, n = i — 1 , 



where we have dropped some irrelevant numerical factors. From Eqs. ( p.3|) , 
(|2.13| ), and (|3.1|) , it follows that when £ is a non-negative integer then the 



Lame Hamiltonian has i + 1 eigenfunctions of the form 

ipix) = cn^x ■ X (^^) , 
Vcn x/ 

where x is a polynomial of degree at most i. Thus we can write 

ip{x) = x{siax,cnx) , (3.2) 

where x is a homogeneous bivariate polynomial of degree i. Similarly, if i 
is strictly positive then there are i additional eigenfunctions of the Lame 
Hamiltonian of the form 

ipi^x) = dnx ■ fj^snXjCnx) , (3.3) 



9 



where r/ is a homogeneous polynomial in two variables of degree Thus, 
when £ is a non- negative integer the Lame Hamiltonian admits 2i+l algebraic 
eigenf unctions of the form (P?^)-(P73|). 

To show that the latter 2i + 1 algebraic eigenfunctions are indeed the 
classical Lame polynomials, we need to study in more detail Eqs. ( |3.2|) -( ^73|) . 
Fortunately, this can be done in a totally systematic way, by applying the 
general theory developed in Ref. for associating a (weakly) orthogonal 



polynomial family to a one- dimensional QES Hamiltonian. According to the 
prescription of the latter reference (cf. [^, Eqs. (30), (33), and (35)), we 
must choose two different roots Zi ^ Z2 of P{z) such that 

(2j-n-l)P'(^i) + 2Q(^i) ^0, Vj>l. (3.4) 

Then the function 

(with C,{x) and /i(C(x)) respectively given by Eqs. ( [^.13| ) and (|3.1|)), is a 
formal solution of the Lame equation with energy E if and only if the set 
{Pj[E)} ..^Q satisfies a certain 3-term recurrence relation. If we normalize 
ipE SO that Pq = 1, the latter recurrence relation implies that {Pj{E)}j^Q 
is a family of (weakly) orthogonal polynomials. Moreover, if the coefficients 
of the recurrence relation satisfy a positivity condition that we shall state 
and verify below, the critical polynomial Pn+i has n + 1 different real 
roots Ei, < i < n. Since it is easily shown that each of these roots is 
also a root of the polynomials Pj{E) with j > n + 1, for i = 0, 1, . . . ,n 
the function ipEi{x) is a genuine solution of the Lame equation ( |1.1| ) with 
energy E = E^. These n + 1 exact solutions can be obtained algebraically, by 
computing the n polynomials Pi{E), . . . , Pn{E) from the recurrence relation, 
and in fact coincide with the n + 1 eigenfunctions of the Lame Hamiltonian 
H algebraically computable by diagonalizing the restriction of H to fiVn (or, 
equivalently, the restriction of H to Vn)- 

It is immediate to verify that ( p.4|) is satisfied for both solutions (|2.19|) 
and ( |2.20| ) simply by taking zi = — i. If we choose Z2 = i, from Eq. ( p.5| ) we 



obtain the following expression (up to a constant factor) for the £ + 1 exact 
eigenfunctions of the Lame Hamiltonian coming from the first algebraization 

dH) 



ipEi{x) = — ^ — Pj{Ei) (cnx + isnx)^~-'(cnx — isna;)-' , < z < 



(3.6) 
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where the i + 1 energies Ei satisfy 

P,+i(E,)=0, 0<i<i. (3.7) 
In the same way, the second algebraization ( p^.20D yields the £ eigenfunctions 



(x) = dnx — ^■j — Pj {Ei){cnx + isnxY ^ -'(cna; — i 



snx 



j=0 



<i< 



(3i 



where i P, > is the orthogonal polynomial family constructed from the 



solution ( p^.20 ), and the £ energies Ei are given by 

Pi{Ei) = 0, 0<i<£-l. 



(3.9) 



We shall next derive the recurrence relation defining the polynomials Pj 
and Pj. According to the general procedure described in Ref. [^], the latter 
relation is determined by the coefficients of the auxiliary polynomials P{w), 
Q{w) obtained by applying to P{z) and Q{z) the projective transformation 



w 



Z — Z\ 
Z- Z2 



(3.10) 



namely 



P{w) 
Q{w) 



{l-wY 
{z,~z,Y 
{l-wf 



p 



Zl - Z2 



Q 



Zl - Z2W 
1 — W 
Zl - Z2W 

1 — w 



(3.11) 
(3.12) 



(with Zl = —i = —Z2 in the present case). From Eqs. ( p.l2| ) and ( 3.11| ) we 
easily obtain the following formula for P, valid for both algebraizations ( |2.19| ) 
and (^) 



P{w) = mw^ — 2 {2 ~ m) w"^ + mw . 



(3.13) 



Similarly, Eqs. (2^) and (|3.12|) imply that the polynomial Q associated to 
the first algebraization ( |2.19| ) is given by 



Q = -m£{w^-l) 



(3.14) 
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while for the second algebraization ( |2.2CI| ) we obtain 

Following Ref. we define new polynomials Pj{E) and Pj{E) by 



(3.15) 



PjiE) 



my (2£-l)!! f PjiE) 
2) (2£-2j-l)!! \p,{E) 



(3.16) 



It is shown in Ref. |23] that the polynomials Pj satisfy the recurrence relation 



= iE- bj) P, - aj Pj_, , J > 



with P_i = and 

a, = ^m'j (2j - 1) (2£- 2j + 1) {£ - j + 1) 

b, = ^mi{i+l) + ^i2-m) (£-2jf . 

Likewise, the polynomials Pj verify the recurrence relation 

P,+i = (E-6,)P,-a,P,_i, j>0, 

where P_i = 0, and the coefficients aj and bj are given by 

1 



a, = -m'j (2j + 1) (2£ - 2j + 1) (i - j) 



' 4 
1 



1 



- ^mi{i + l) + -i2-m)ii- 2j - ly 



(3.17) 

(3.18) 
(3.19) 

(3.20) 

(3.21) 
(3.22) 



Note that from the normalization Pq = Pq = 1 and Eq. (|3.16|) it follows that 

Po=h = l- (3.23) 

Furthermore, Eqs. ( p.l8| ) and ( p.21| ) imply that the coefficients aj and cij 
are strictly positive for 1 < j < £ and 1 < j < ^ — 1, respectively. By 
Lemma 1 (Sec. 1.8) of Ref. 0], this positivity condition guarantees that all 

the roots of the critical polynomials Pe+i and Pe are real and simple, as we 
had anticipated. 

Our next task is to simplify Eqs. ( ^.6| ) and ( |3.8| ) using the properties of 
the polynomials Pj and Pj. We shall start by proving the following useful 
result: 
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Proposition 3.1. The algebraic eigenfunctions ( p.6[ ) and ( |3.8| ) are always 
either real or purely imaginary. 



Proof. Let us rewrite Eqs. ( |3.ti| ) and ( |3.^ ) in terms of the elliptic amplitude 
ip{x) = amx, defined by 

cnx + i snx = e^'^'-^'' . 

We thus obtain the expressions 
i 

^eA^) = J]p.,e'(^-2^>(^), 0<t<e, (3.24) 

and 

e-i 

(x) = dn X ^ prs ei(^-2^-i)^W , < r < £ - 1 , (3.25) 

s=0 

where we have set, for convenience, 

Prj = ^-y^ PJ{E^), Prs = Ps{Er) • (3.26) 

Let us concentrate, for definiteness, on the eigenfunctions of type (|3.24| ). We 
rewrite Eq. (|3.24|) in the form 

^eA^)= E N e'^'"''^^^'^ + PM-i (3.27) 

0<j<i/2 

where we take Pi//2 = when i is odd. Taking into account that, by the 
recurrence relation ( |3.17|) - (|3.19|) , all the coefficients real, we conclude 

from the Eq. {^727\) that 



Re V'i?,(a;) = ^ {pij + Pi/^j) cos(£ - 2j)ip{x) + p,-^^ , (3.28) 

0<j<i/2 

Iui^eA^) = {Pij -Pi/-j) sm{i-2j)(.p{x) . (3.29) 

0<j<i/2 

Let us now introduce the polynomial family 

n,iE) = ^P,(i?), (3.30) 



13 



in terms of which the coefficients pij are simply expressed by 

Pij = TTjiEi) , <i,j <£. 
Defining two additional families of univariate polynomials <yj{E) and Pj{E) 

by 

a,{E) = 7^,{E) + 7r,„,(E) , p,{E) = n,{E) - 7r,_,(E) , (3.31) 

from Eqs. we have 

ReV^ij, = ^ <yj{Ei) = () for < j < [£/2] , 
Im^E^ = pj{Ei) = for < j < 1/2 

(where [ ■ ] denotes the integer part). From Eqs. ( p.l6| )- (|3.19|) and ( p.30|) , 
it immediately follows that the polynomials T^j{E) satisfy the three-term 
recurrence relation 

(j + 1)(2£ - 2j - l)7r,+i = ^{h, - E)7T^ - (2j - 1)(^ - j + l)7r,_i (3.32) 

with the initial conditions 

7r_i(E)=0, 7ro(E) = l, (3.33) 

where the coefficients bj are defined by Eq. (|3.19|) . It is easy to see that 
Eq. (|3.31| ) implies that the polynomial families {crj^E)}^^^ and {pj{E)}j^^ 
satisfy the same three-term recurrence relation (|3.32|) as the family {TTj{E)}j^Q. 
Moreover, from Eq. (|3.7|) we obtain 



a_i{Ei) = p_,{Ei) = , 0<i<i. (3.34) 

From this, and the fact that the coefficients {crj{Ei)} .^^ and {pj{Ei)} 
satisfy a three-term recursion relation of the type ( ^.3^1 ), we conclude that 
the vanishing of (Jo{Ei) or of po{Ei) automatically implies the vanishing of 
aj{Ei) or of pj{Ei), respectively, for all values of j > 0. Thus we can write 

Re^JE, = ^ ao{Ei) = , Im^^^;^ = ^ po(^i) = . (3.35) 

To complete the proof of the Proposition, we simply note that Eqs. (|3.33|) - 
( |3.34| ) and the fact that {aj{Ei)}j^Q and {pj{Ei)}j^^ satisfy the same three- 
term recursion relation as {TTj{E)}jyQ imply that 

a,{Ei) = ao{Ei)7r^{Ei), j>0. 



14 



and, in particular 



^ ao{Ei) [1 - TTe{E,)] = a^{Ei) poiE^) = . (3.36) 

The proof for the eigenfunctions of type ( |3.25| ) is totally analogous, and will 
therefore be omitted. □ 

Let us order, from now on, the i+1 algebraic eigenvalues Ei of the Lame 
equation in increasing order, i.e., 

Eq < El < ■ ■ ■ < Ei_i < Ei , 

and similarly for the i eigenvalues Ei. From Lemma 1 (Sec. 1.8) of Ref. |^ 
and the positivity of aj for 1 < j < £ we easily deduce that 



sign [{-lyPeiE,)] = i-iy 



(3.37) 



which implies (by Eq. ( |3.30|) ) that vr^(-Ei) has the sign of (—1)*. Thus for 
even i we have a^^Ei) = 1 + 7ii{Ei) > 0, and from ( p.36| ) we conclude that 
Po{Ei) = 0. In the same way we establish that cTo{Ei) = for odd i. From 
Eq. ( p. 351 ) we obtain 



Reip 



E2r 



Im if) 



E2r 



0. 



(3.38) 



We can therefore write the £ + 1 exact eigenfunctions of type ( |3.(j| ) of the 
Lame equation as the Jacobi-Fourier series 

(j)2r{x) = ^E2r{^) = ^2r,i COs(£ - 2j)ip{x) + Psr.f (3.39) 
0<i<£/2 

02s+i(a;) = ^E2s+iix) = Yl sin(£ - 2j)^(x) , (3.40) 



0<j<i/2 



with 



C"2rJ = 0'j{E2r) — P2r,j + P2r,l-j , P2s+l,j = Pj{E2s+l) — P2s+l,j — P2s+l,l-j , 

(3.41) 



and 



< r < 



< s < 



- 1" 



(3.42) 
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We can derive in a totally analogous way an equivalent formula for the ^ 
eigenf unctions of type ( p.27| ), namely 



2r{x) = ^E^M) = dnx 



0<i<|(^-l) 



<p2s+i{x) = ipE2s+i(^) = ^ /^2s+ij sin(£ - 2j - l)ip{x) 

o<i<|(^-i) 

where now 

d'2r,j = P2r,j + P2r,i-j-l , P2s+l,j = p2s+l,j — p2s+l,l-j-l , 

and 



(3.43) 
(3.44) 



< r < 



< s < 



- 1 



(3.45) 



(3.46) 



Our next step is to further simplify Eqs. ( |3.39|) -( p.46|) for the algebraic 
eigenfunctions of the Lame Hamiltonian. To this end, recall that by con- 
struction the Chebyshev polynomials Tj and Uj (j = 0, 1, 2, ... ) satisfy the 
identities 



cos(ja) = Tj(cosQ;) , 
sin(ja) = sin a Uj^i{cos a) 



(3.47) 
(3.48) 



where the latter equality is formally valid for j = if we set U-i = 0. It 
follows that 



(p2rix) = X a2r,j r£_2j(cnx) +P2r,| 
0<j<i/2 

Wi(a;) = snx X p2s+i,jUe^2j-i{cnx) 
o<j<e/2 



(3.49) 
(3.50) 



and, similarly. 



■>2r 



'x) = diax P2r,i(<?-i) + XI cr2r,jTi_2j -i{cnx) 

o<i<|(f-i) 

l2s+i{x) = snx dnx ^ P2s+i,j Ui-2j-2{cnx) . 

0<j<^(i-l) 



(3.51) 
(3.52) 
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It is now straightforward to check that Eqs. ( |3.49| )-( |3.52| ) encompass the 
eight types of classical Lame polynomials listed in Ref . @] . Indeed, note first 
of all that the polynomials Tj{t) and Uj{t) have the parity of (—1)-' under 
the reflection x i— > —x. Therefore the formulas 

T,,(t) = ^^, t/,,(t) = ^^; J =0,1,2,.... (3.53) 

define even polynomials T2j{t) and U2j{t) of degree 2j. If £ = 2N {N = 
1,2,...), Eqs. ( |3.49| )-( p.52| ) yield the following four types of polynomial — 
in sn, cn, dn — solutions of the Lame equation: 

N-l 

(f)2rix) = P2r,N + ^ Cr2rJ T2Af-2j(cn x) = ttrisn^ x) (3.54) 
j=0 
N-l 

<p2s+i{x) = snxcnx p2s+i,j f^2Af-2j-2(cnx) = snxcnx/5s(sn^x) (3.55) 
i=o 

Af-l 

02s (a;) = cnxdnx a2sj T2Ar_2i-2(cnx) = cnxdnxas(sn^x) (3.56) 

j=0 
N-l 

(t>2s+i{x) = snxdnx'^ p2s+i,jU2N-2j~2{cnx) = sna;dnx/3s(sn^x) , (3.57) 
i=o 

where 

0<r<A^, 0<s<A^-l, (3.58) 
and arit), (3 sit), cxsit), and 13 sit) are polynomials of respective degrees 

deg ar = N , deg (3s = deg ois = deg I3s = N — 1 . (3.59) 
Similarly, if £ = 2N + 1 (A^ = 0, 1, 2, ... ) the exact solutions of the Lame 
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equation ( |3.49| )-( |3.52| ) can be written as follows: 



N 

(p2r{x) = cn X (72r,j ^2jV-2j {cn x) = cn x 7,, (sn^ x) (3.60) 

j=0 

N 

^2^+1 (a;) = sn X p2r+i,i t^2Ar-2i(cnx) = snx5r.(sn^x) (3.61) 
i=o 

N-l 

J2rix) = dna; p2r,N + ^ ^2r,j ^27v-2j(cnx) = dnx7r(sn^a;) (3.62) 

3=0 
N-l 

2 



^2s+i(a:) = snxcnxdnx p2s+ij t/2Ar-2i-2(cnx) = sna;cna;dnx5s(sn^ 

j=0 



X) 



(3.63) 



where the range of the indices r, s is still given by Eq. ( |3.58| ), and 7r(t), Sr{t), 
7r.(t), and 6s{t) are polynomials with 



deg 7r = deg 6r = deg •jr = N , deg 6s = N — 1 . 



(3.64) 



3.1 Examples 

The explicit expressions for the Lame polynomials in terms of Chebyshev 
polynomials derived at the end of the previous section can be used to algo- 
rithmically compute the Lame polynomials and their corresponding energies 
for low values of i. 

When i = I, Eqs. ( |3.6(]| )- (|3.63|) for = immediately yield the following 
well-known formulas for the Lame polynomials of order 1: 



(x) oc ciax, 



(a;)ocsna;, ^oocdna;. 



(3.65) 



The critical polynomials P2 and Pi are easily computed using the recurrence 
relations (^J7\)-{^^, with the result 



p2{E) = - {m + 2)E + m + l, Pi{E) = E - m . 



(3.66) 



The energies associated to the eigenfunctions (|3.65|) are just the roots of each 
of the latter polynomials ordered increasingly, namely 



El = m + 1 ] Eq = m . 



(3.67) 
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For i = 1 the critical polynomials are 

Py,{E) = E'^ - (5m + 8)E^ + ^{m^ + 8m + 4)E - 12m (m + 4) , (3.68) 

P^{E) = E"^ - (5m + 2)E + (m + l)(4m + 1) . (3.69) 

The energies corresponding to the Lame polynomials of order two are there- 
fore given by 



E- 



°| = 2^1 + m=F Vm? — m + ij , Ei = A + m; 



(3.70) 



Eo = m + l, Ei = 4m + 1. (3.71) 



Note that the factor of 2 in (131701 ) is missing in Ref. [Q, p. 205. The Lame 

polynomials are easily computed using (|3.54|) - (|3.57|) with = 1, with the 
following result: 

= sn^ X (l + m± Vm^ — m + 1^ , (3.72) 

3m V / 

^ = snxcna;; (3.73) 

^ = cnxdna;, ^ = snxdna;. (3.74) 

When i = 3 the critical polynomials are given by 

Pi{E) = [E^ -2(5 + 2 m)E + 3(3 + 8 m)] 

X [E^ - 10(1 + m)E + 3(3 m^ + 26 m + 3)] , (3.75) 

PsiE) = [E- 4(m + 1)] [E^ - 2(2 + 5m)E + 3m(3m + 8] , (3.76) 

from which the following exact energies are computed: 



^0 

E2 

El 
E3 

Eo 
E2 



2m + 5 + 2 Vm2 - m + 4 , (3.77) 



5(m + 1) + 2 V4m2 - 7m + 4 ; (3.78) 



= 5m + 2 + 2V4m2 - m + 1 , (3.79) 
El = 4(m + 1) . (3.80) 
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From Eqs. ( |3.6CI| )-( |3.63| ) with = 1, the corresponding Lame polynomials 
are given by 



cnx 



sn^ X 



(m + 2 ± V ra^ — m -\- A] 

5m V / 



(3.8i: 



snx 



sn^ X 



— (2{m + 1) ± V4m2 - 7m + 4 
5m V 



(3.82) 



dnx 



sn^ X 



— f 2m + 1 ± 74 
5m V 



m^ — 171+ 1 



(3.83) 



— = sn X cn X dn x . 



(3.84) 



The latter formulas for the Lame polynomials of order three and their corre- 
sponding energies are seen to coincide with those given in Ref. 0]. 

When i is greater than three, it becomes increasingly more difficult to 
compute the Lame polynomials of order £ and their corresponding eigen- 
functions in closed form. This is essentially due to the fact that the high 

degree of the critical polynomials Pi and P^+i in general makes it impossible 
to exactly compute the energies Ei and Ei, which enter in the definition of 
the coefficients pij and pij. For example, for i = 4 the critical polynomials 

P5 andP4 are given by 

P^^E) = [E^ - 10 E (2 + m) + 64 + 136 m + 9 m^] 

X [p3-20E2(l + m) + 16E(4 + 21m + 4m2)- 640m(l + m)] , 

(3.85) 

Pi{E) = [E^ - 10 E (1 + m) + 9 + 46 m + 9 m^] 

X [E^ - 10 E (1 + 2 m) + 9 + 136 m + 64 m^] . (3.86) 

It is not difficult to verify that the roots of the quadratic factor of P5 are the 
energies Ei and E-^, so that we have the exact formula 



E->, 



5(m + 2) =F V4m2 - 9m + 9 . 



(3.87) 



20 



The corresponding Lame polynomials are easily computed using Eq. ( p.55| ): 



snx cnx 



sn^ X — - — f 2m + 3 ± V Ara^ — 9m + 9 
7m \ 



(3. 



Similarly, the fact that P4 factorizes into two quadratic polynomials allows 
us to exactly compute the energies Ei in this case, obtaining 



Eq 

E2 
El 
E3 



5(m+ 1) T2V4m2 + m + 4, (3.89) 

5(2m +1)t 2\/9m2 -9m + 4. (3.90) 
By Eqs. ( p.56| )-( ^.57| ), the corresponding eigenfunctions are given by 



cnx dnx 



sn^ X 



(2(m + 1) ± V4m2 + m + i] 
7m V / 



(3.91) 
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snx dnx 



sn^ X 



f 5m + 8 ± AV9m^ - 9m + 4 
28 m V 



. (3.92) 



The remaining Lame polynomials of order four can be expressed using Eq. ( p. 54]) 
in terms of the three roots Ei [i = 0, 2, 4) of the cubic factor of P5, with the 
following result: 



— = sn X 16(m + 1) — EA sn x 

16 14 m 



+ 



1 



280 m^ 



[E^ - 20(m + l)Ei + 6Am^ + 296m + 64] . (3.93) 



For instance, for m = 1/2 the roots £'0, E2 and £'4 can be computed in closed 
form: 



E, 



iot'^VTs, ^2 = 10. 



(3.94) 
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Substituting into Eq. ( |3.93D we obtain the following exact expressions for the 
corresponding Lame polynomials for m = 1/2: 



sn^x-'^(7± Vltj sn^ X + ^ (^4 ± v^) , (3.95) 



— = sn^ X — 2 sn^ x -\ — . (3.96) 
16 5 

4 Case II: ^ is a positive half-integer 

This is the most interesting case, since we shall see that it leads to the non- 
meromorphic Lame functions, whose algebraization has proved more difficult 
than that of the Lame polynomials. 

The gauge factor can be computed from Eqs. ( p.9|) , ( p.l3|) , and (|2.21|) - 
(12:2^ ), with the result 



fi {({x)) = cn" x{cnx±ik' sn x) ' . (4.1) 

The upper and lower signs in this formula correspond, respectively, to the 
solutions (12:211) and j^^ , 



n = i-l, (4.2) 



and we are using the customary notation 



k = \/m, k! = \/\ — m . (4.3) 

Thus when £ is a positive half-integer the Lame equation has 2n + 2 = 2^ + 1 
algebraic solutions of the form 

?/'=^(a;) = cn"a; (cnx ± i/c'snx)^/^ j , (4-4) 

where is a polynomial of degree at most n. 

We shall now use the same strategy as in the previous section to simplify 
the eigenfunctions ( [4.41 ). It is straightforward to check that in this case only 
the two roots ±i/A;' of the polynomial Pi^z) in (|2.12|) satisfy ( p.4|) . For reasons 
that shall become clear later, we shall choose 

i i 
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where again the upper sign corresponds to ( |2.21| ) and the lower sign to ( p.22[ ). 
Using the same convention for the eigenfunctions ^p^± and the polynomials 

i 

Pp, a straightforward calculation yields 
i)^±{x) = ^^-jLp±(^E^) {cnx±ik'snx)''-^+^cnxTik'snxy . (4.5) 

We have dropped an inessential constant factor in the latter formula, and 
the energies Ef are defined by 

Pn+iiEt) = 0, 0<z<n. (4.6) 

(We shall prove below that the roots of -P^i are indeed real and simple.) 

We now compute the polynomials P^, or more precisely the recurrence 
relation defining them. In principle, we should use the notation P^{w) and 
Q^(tf) to denote the polynomials P{w), Q{w) in Eqs. ( |3.11| )-( |3.12| ) corre- 
sponding to each of the two cases ( p.21j ) and ( p.22|) . Remarkably, a direct 
calculation shows that this is actually unnecessary, since 

P+{w) = P'{w) = P{w) = - [mw^ + 2 {2 - m) w'^ + mw] . (4.7) 

-s -s 777, X 

Q+{w) = Q'{w) = Q{w) = —{n + 2)w^ + {2-m)w - -mn . (4.8) 

This implies that the recurrence relations determining the two sets of or- 
thogonal polynomials \^P^{E)^ are in fact the sameQ. Since we have 

normalized both sets so that P^{E) = 1, it follows that both polynomial 
families are identical, namely 

P^iE) = Pr(E) ^ P^iE), J = 0, 1, 2, ... . (4.9) 

In particular, the above equality for j = n + 1 and Eq. ( [4. 61 ) show that 

Et = Er =Ei, 0<i<n. (4.10) 

Moreover, to each of these n + 1 energies correspond two algebraic eigenfunc- 
tions of the form (|4.5|), which on account of ([4.9|) can be written as 



" (—iV 

= J2^-j^PjiEi) (cna;±iA;'snx)"-^+5(cnxTiA;'snx)^ (4.11) 



■^From Eq. (2.1J) we have c+ = , and from the fact that R is invariant under projec- 
tive transformations, [|l9|, and Eqs. ( p.6D and ( [4. 7] ) we obtain 

ct = -^{n + 2) (coo - c+ ) + c+ = c" . 
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From Eqs. ( [4.7D -( ^y8D we see that the coefficients of the polynomials P{w) and 
Q{w) are all real. This implies, that the coefficients of the recurrence 
relation determining the polynomials Pj{E) are also real. Since Po{E) = 
1, it follows that Pj{E) is a polynomial with real coefficients, for all j = 
0,1,.... From Eq. ( [4.11|) and the positivity of the coefficient aj for 1 < j < n 
(cf. Eq. ([4.14|) below), it follows that all the energies Ei {0 < i < n) are real 
and simple, and therefore 

'^E,=i^t^ 0<i<n, 
where the overbar denotes complex conjugation. If we set 

ipi{x) =ip^^{x) , <i <n, 

then for each root Ei of the Lame equation has two algebraic rea/-valued 
solutions given by 

(j)l{x) = Re^,{x) , (f)^,{x) = Im^iix) . (4.12) 

The recurrence relation satisfied by the polynomials Pj{E) can be easily 
written down from (|4.7| ) - (|4.8| ) using the prescription of Ref. [^]. Indeed, if, 
following the latter reference, we write 



j! \kj (2j)! • 

then the polynomials Pj satisfy a three-term recurrence relation of the form 
(|3T^ ), with 

a, = ^ m^j{2j - l)(2n - 2j + 3)(n - j + 1) , (4.14) 



bj = 1 (2n + l)(m + 2ra + 1) - (2 - m)j {2n - 2j + 1) . (4.15) 



We are ready to show that, when £ is a half- integer, the 2i + 1 algebraic 
solutions of the Lame equation of the form ( [4.12| )- (|4.11|) obtained in the 
previous section are precisely the non-meromorphic Lame functions. 

To this end, we need to simplify Eq. (|4.11| ). In the ffist place, since 

|cnx ± i/c'snx| = dnx (4.16) 
(recall that dnx is positive for all real x), we can write 

cnx ± i fc'snx = dnxe^^'^^"") , (4.17) 
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so that 



cos Qix) = , sin e(x) = k' . 4.18 

dn X dn X 



Secondly, we have 

1/2 ^ _^ 

'V2 



cnx + i fc' sn x)^^^ = ^Vdnx + cnx + ie(x)Vdnx — cnxj , (4.19) 



where we have set 

e(x) = sign(snx) . (4.20) 

Note that the functions Vdn x + cn x and e(x) Vdn x — cn x are C°° on the 
whole real line except at the points x = 2(2r + 1) K. We shall also find useful 
in the sequel the identities 

r- 1 dn X ± cn X r- 

vdnx ± cnx = e(x) — ; vdnx^pcnx. (4-21) 

k' snx 

From ( [4.17|) it follows that 

n 

"^^Pij (cnx + i k' snx)"~-'(cnx — i k' snx)-* 

j=0 

n 

= dn" X ^ Pij e' , (4.22) 

j=o 

where the coefficients pij are again defined by Eq. ( p.26|) . The sum in the 
latter equation is of the same time as ( p.24| ), and therefore can be expressed 
as follows: 

n 

i=0 0<j<n/2 

+ isin6'(x) ^ pij Un^2j -i [cos 9 (x)) 

0<i<n/2 

= Pi,n/2 + 5Z ^) 
0<j<n/2 

+ iA;'sdx [/„_2j-i(cdx) (4.23) 

0<j<n/2 

where we have used the customary abbreviations 

cnx snx 

cd X = , sd X = , 

dn X dn X 
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and we have again set 



Pij ~l~ Pi,n—j 



Pij Pij Pi,n—j ■ 



(4.24) 



From Eqs. (|4.19|) and (^4.21|) (dropping the constant factor ±l/-\/2) we easily 
obtain 



'>l{x) = Vdn X + cn a; dn^a; 



Pi,n/2+ ^ CrijTn-2j{cd 
Q<j<n/2 



X 



+ (cdz - 1) pij f/„_2j_i(cd; 

0<jr<n/2 



. (4.25) 



An analogous calculation leads to the following expression for the second 
solution (ti^ix): 



4)1{x) = e{x) Vdn x — cn x dn"x 



Pi,n/2 + 

Q<j<n/2 



(TijTn-2j{cdx) 



+ (cdx + 1) ^ Pij f4_2j_i(cdx) 

0<i<n/2 



. (4.26) 



Let us now study in more detail the properties of the 2i+l solutions of the 
Lame equation of the form ([4.25|) -( fi.26|) that we have just obtained. From the 
differentiability properties of the functions Vdn x + cn x and e(x)Vdnx — cnx 
mentioned above, it follows that 4>}{x) and (pf{x) are C°° solutions of the 
Lame equation (|1 . 1| ) with E = Ei in the open interval (— 2^, 2^), and are 
clearly continuous on [—2K,2K]. The product of dn"x with the functions 
within square brackets in Eqs. ( [4.25[ ) and ( |4.26| ) defines two C°° functions 
on the whole real line (they are in fact polynomials in snx, cnx, dnx), with 
(real) period 4K. On the other hand, it is immediate to check that if we pro- 



long /(x) = -\/dn X + cn X or /(x) = e(x) Vdnx 
[—2K, 2K] anti-periodically, i.e. requiring that 

/(x + 4if) = -/(x). 



cnx outside the interval 



(4.27) 



then / becomes a C°° function. (Indeed, in a neighborhood of ±2K the 
prolongations of Vdn x + cn x and e(x) Vdnx — cnx are equal, respectively, 
to ±k' snx/Vdnx — cnx and ±a/ dnx — cnx, the latter functions being ob- 
viously C°° at ±2K.) It follows that, if we define ^^(x) and ^^(x) outside 
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\—2K,2K] by the prescription ( |4.27| ), then and (f)j are C°° on the whole 
real line, and have (real) period 8 K. Furthermore, from Eqs. ( [4.25|) and 
( |4.2(j| ) it follows that (fl^x) and 0^(a;) are, respectively, even and odd func- 
tions of the variable x. It is also clear from the well-known equalities 

sn(x + 2K) = — sn a; , cn{x + 2K) = — cnx , dn{x + 2K) = dn a; , 

that (pl and (pj are related by 

<P^{x) = ^l{x + 2K). (4.28) 

From this equality and the anti-periodicity condition ([4.27|) it follows that 

[<Plix + 2K) ±icpj{x + 2K)\ = Ti [<PKx) ± i0^(x)] , (4.29) 

so that the quasi- momentum associated to each of the 2i+l algebraic energies 
Ei is equal to n/^AK). 

When n is odd, it is convenient to rewrite Eq. ( 4.25| ) in the following way: 

|(n-l) 

Cdx ^ (7ijfn-2j~l{cdx) 



(f)l{x) = Vdn X + cn X dn'' 



.X) 



Vdn a; -|- cn a; dn" 



j=0 

i(n-l) 

+ {cdx -1) ^ pij f/„_2j_i(cd: 

j=0 
i{n-l) 

cna; ^ (Jij T„_2j_i(cdx) 

j=0 

|(n-l) 

+ (cn a; - dn x) ^ Pij f/„-2j-i (cd x) 

3=0 



X 



(4.30) 



Taking into account the parity properties of the Chebyshev polynomials, and 
the fact that cn^ x and dn^ x can be expressed linearly in terms of sn^ x, it 
follows that 



0,-(x) = Vdnx + cnx [cn x ai (sn^ x) + dnx Pi{sn^ x)~\ , (4-31) 
where ai{t) and l3i{t) are polynomials in t of degree |(n — 1) given by 

|(n-l) 

ai(sn^a;) = dn''"^^ ^ jdij r„„2i-i(cdx) + [/„_2j-i(cdx)| , (4.32) 

j=0 
|(n-l) 

^(sn^a;) = -dn"-^a; p,,- f/„_2,_i(cda;) . (4.33) 

i=o 
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Comparing ( ^.311 ) with Ref. (p. 97), we see that cp} is proportional to the 



ri + - 

non-meromorphic Lame function Ec^' ^ , 



0<i<n=i- - 
- - 2 



(4.34) 



where rj is the number of zeros of (J)} in the interval (0, 2K). From Eq. ( ^4.26| 
or, alternatively, from Eqs. ( [4.28|) ) and ( [4.31|) , we also obtain 

0^(a;) = e{x)Vdi]lc'^^cnx [cnxaj(sn^x) — dnx (3i{sn^ 

so that 

1 



X 



h ix) oc Es 



I 1 



< i < n 



(4.35) 
(4.36) 



Similarly, when n is even we rewrite Eq. (^.25|) as follows: 



(j)\{x) = Vdn X + cn X dn"x 



Pi,n/2 + ^ <7ij Tn-2j{cd 



XI 



j=0 



Cdx (cdx - 1) ^ Pij Un-2j-2icd 



X] 



j=0 



Vdn X + cn X dn" ^ 



X 



2 

dn^ X (^Pi,n/2 + ^ O-ij Tn-2j{cd 



j=0 

-1 



+ cn X (cn X — dn x) pij f/n-2j-2(cd 

j=0 



X 



(4.37) 



We thus obtain 



(plix) = Vdn X + cn X [Q;j(sn^x) +cnx dnx/3j(sn^x)] , (4.38) 

where ai(t) and Pi{t) are polynomials in t of respective degrees | and f — 1, 
defined by 

Pi,n/2 + ^ o-jj r„_2j(cdr 



aj(sn^ ^) = dn" x 



[X 



j=Q 



--1 
2 



A(sn^ x) 



+ dn" ^ X cn^ x p^j f/„-2i-2(cdx) , 

j=0 

\(n-l) 

-dn"~^X ^ PijUn-2j-2{cd: 
j=0 



Ix 



(4.39) 
(4.40) 
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From ( ^.28| ) we have 

= e(x) Vdn X — cn X [Q;i(sn^x) — cnx dnx/3j(sn^ 



X 



(4.41) 



Comparing, again, with Ref. we deduce that Eqs. ( [4.34|) and ([4.36|) also 
hold when n is an even non-negative integer. 



4.1 Examples 

The case n = 0, i.e, I = 1/2, is particularly easy to deal with in the framework 
of our formalism. Indeed, in this case we can exactly compute one eigenvalue 
Eq of the Lame Hamiltonian, its corresponding two linearly independent 
eigenf unctions given by (f4.38| )-( ^.41|) with n = 0: 

0q(x) = Vdn X + cn X , 0o(^) ~ e(x)Vdnx — cnx . (4.42) 

The eigenvalue Eq is the root of the first degree polynomial Pi or, equiva- 
lently, of Pi (since both polynomials differ by a constant factor, by (|4.13|) ). 
From Eq. (|3.171 ) with j = we get 

Pi(^) = (^-M ^0 = ^-^0, 
and from ( [4.15|) (with n = 0) it follows that 

Eo = ho = \{l + m) = ^{l + e). (4.43) 

We have thus shown in a purely algebraic fashion that the general solution 
of the Lame equation 

^"(x) + ^ (l + A;^ - 3 A;^ sn^ x) ^(x) = 

is a linear combination of the two functions ([4.42|) . 

Similarly, for n = 1, i.e, £ = 3/2, the two algebraic eigenvalues of the 
Lame Hamiltonian are the roots of the polynomial 

P2(E) = - - (m + 1) E + — (3 + 23 m + 3) , (4.44) 
2 16 

namely. 



^ K = -{m + l)T -m + l. (4.45) 
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Their associated eigenfunctions are easily computed using Eqs. ( [4.31| )-( P:.33D , 
which for n = 1 reduce to 



-P,{E,)-l = --P^{E,)-l = - 
m m 



-{m + 3)-Ei 



Ui = 2pio = 2 
A = Pii - Pio 

From Eq. ( |4.31| ) we obtain 

= Vdn X + cn X mcnx + (1 — m±V — m + 1 ) dn x 



1. 



(4.46) 



where we have omitted a trivial constant factor. The remaining two eigen- 
functions 0Q and (pl can be easily obtained from (|4.46| ) using ( [4.35| ): 



ex 



) 



X — cnx 



m cnx 



1 — m ± V — m + 1 ) dn x 



(4.47) 

The latter formulas are easily seen to coincide with the ones at the end of 
Ref. §. 

As in the case of integer i, the computation of the non-meromorphic Lame 
functions in closed form quickly becomes unmanageable as the degree of the 
critical polynomial Pn+i increases. For example, for i = 5/2 (i.e., n = 2) the 
algebraic eigenfunctions of the form ( |4.371 ) can be expressed as follows 

4>} 



24 m2 



Vdnx + cnx { -25 - 130 m + 87 + (104 + 40 m) Ei - 16 E^ 



+ m sn^ X [25 + 34 m + 9 - (104 + 40 m) Ei + 16 Ef] 

+ 2 cnx dnx [25 + 430 m + 21 - 8 (13 + 11 m) + 16 Ef] } , 

(4.48) 

where E^ is one of the three real roots of the critical polynomial 



35 



7 



Ps(E) = E'^ - ^{m + 1)E^ + — (37m'^ + 138m + 37)^ 
4 16 



64 



m 



+ 1) (45m2 + 794m + 45) . (4.49) 

The remaining three eigenfunctions (pf {i = 0,1,2) are obtained from the 
Eq. (|4.48| ) using Eq. (|4.41| ). When m = 1/2, the energies Ei can again be 
computed explicitly, i.e.. 



Eq 
E2 



35 a; 

irT V7, 



(4.50) 
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Consequently, Eq. ( [4.4^ ) for the eigenfunctions yields the following explicit 
expressions: 



3.1 




Vdn X + cn X (5 ± V7) sn^ x — 2 (2 ± V?) cn x dn x — 7 =F 2V7 



1^1 
16 



= Vdn X + cn X ^sn^ x + 2 cn x dn x — . 



(4.51) 
(4.52) 



References 

[1] Ince E L 1939-40 Proc. Roy. Soc. Edmb. 60 47-63 

[2] Ince E L 1939-40 Proc. Roy. Soc. Edmb. 60 83-99 

[3] Erdelyi A, Magnus W, Oberhettinger F and Tricomi F G 1955 Higher 
Transcendental Functions, vol. Ill (New York: Mc-Graw Hill) 

[4] Arscott F M 1964 Periodic Differential Equations, (Oxford: Pergamon) 

[5] Alhassid Y, Giirsey F and lachello F 1983 Phys. Rev. Lett. 50 873-6 

[6] Braibant S and Brihaye Y 1993 J. Math. Phys. 34 2107-14 

[7] Kofman L, Linde A and Starobinsky A A 1994 Phys. Rev. Lett. 73 
3195-8 

[8] Greene P B, Kofman L, Linde A and Starobinsky A A 1997 Phys. Rev. 
D 56 6175-92 

[9] Turbiner A V 1989 J. Phys. A: Math. Gen. 22 L1-L3 

[10] Gonzalez-Lopez A, Kamran N and Olver P J 1994 Contemporary Math- 
ematics 160 113-40 

[11] Ushveridze A G 1994 Quasi-Exactly Solvable Models in Quantum Me- 
chanics (Bristol: lOP) 

[12] Olver P J 1997 GR0UP21: Physical Applications and Mathematical As- 
pects of Geometry, Groups, and Algebras., ed H-D Doebner et al. (Sin- 
gapore: World Scientific) 285-95 

[13] Razavy M 1980 Am. J. Phys. 48 285-8 

31 



[14] Finkel F, Gonzalez-Lopez A, and Rodriguez M A 1999 Preprint [math- 
I ph/9905020| 

[15] Ward R S 1987 J. Phys. A: Math. Gen. 20 2679-83 

[16] Ulyanov V V and Zaslavskii O B 1992 Phys. Rep. 216 180-251 

[17] Brihaye Y and Godart M 1993 J. Math. Phys. 34 5283-91 

[18] Turbiner A V 1988 Commun. Math. Phys. 118 467-74 

[19] Gonzalez-Lopez A, Kamran N and Olver P J 1993 Commun. Math. 
Phys. 153 117-46 

[20] Turbiner A V 1992 J. Phys. A: Math. Gen. 25 L1087-L93 

[21] Finkel F and Kamran N 1998 Adv. Appl. Math. 20 300-22 

[22] Turbiner A V and Ushveridze A G 1987 Phys. Lett. A126 181-3 

[23] Finkel F, Gonzalez-Lopez A and Rodriguez M A 1996 J. Math. Phys. 
37 3954-72 



32 



